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We study the gravitational collapse of compact objects in the Brane- World. We begin by arguing 
that the regularity of the five-dimensional geodesies does not allow the energy-momentum tensor 
of matter on the brane to have (step-like) discontinuities, which are instead admitted in the four- 
dimensional General Relativistic case, and compact sources must therefore have an atmosphere. 
Under the simplifying assumption that matter is a spherically symmetric cloud of dust without 
dissipation, we can find the conditions for which the collapsing star generically "evaporates" and 
approaches the Hawking behavior as the (apparent) horizon is being formed. Subsequently, the 
apparent horizon evolves into the atmosphere and the back-reaction on the brane metric reduces 
the evaporation, which continues until the effective energy of the star vanishes. This occurs at a 
finite radius, and the star afterwards re-expands and "anti-evaporates" . We clarify that the Israel 
junction conditions across the brane (holographically related to the matter trace anomaly) and the 
projection of the Weyl tensor on the brane (holographically interpreted as the quantum back-reaction 
on the brane metric) contribute to the total energy as, respectively, an "anti-evaporation" and an 
"evaporation" term. Concluding, we comment on the possible effects of dissipation and obtain a 
new stringent bound for the brane tension. 
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I. INTRODUCTION 

It is well known that black holes are unstable in four 
(and higher) dimensions because of the Hawking ef- 
fect that is the quantum mechanical production of 
particles in strong inhomogeneous gravitational fields. It 
is also well known that such an effect leads (and is deeply 
linked) to the trace anomaly of the radiation field on the 
black hole background 0, However, for a complete 
description, the semiclassical Einstein equations should 
be solved including the back-reaction of the evaporation 
flux on the metric, which turns out to be an extremely 
hard task (for a recent attempt to incorporate the effect 
of the trace anomaly see Ref. 0). 

In the context of the Randall-Sundrum (RS) Brane- 
World (BW) models @, it was shown in Ref. 6] (see also 
Ref. Q for some recent generalizations) that the collapse 
of a homogeneous star leads to a non-static exterior, con- 
trary to what happens in four-dimensional General Rela- 
tivity (GR), and a possible exterior was later found which 
is radiative Q . If one regards black holes as the natural 
end state of the collapse, one may conclude that classical 
black holes in the BW should suffer of the same problem 
as semiclassical black holes in GR: no static configura- 
tion for their exterior might be allowed. 
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In particular, it was shown in Ref [9j, that all known 
black hole-like metrics on the brane lead to Weyl anoma- 
lies with a natural interpretation in the context of the 
holographic analogy [l(J]- Moreover, such anomalies 
could be related with an instability, as those metrics 
do not seem to have the correct weak field exp ansion 
in RS (for a discussion of this issue see Ref. [llj ) . Fur- 
ther, forcing a static exterior, a trace anomaly outside a 
homogenous and isotropic collapsing star appears which 
is of the same form, but with opposite sign, as that of 
semiclassical black holes. This suggested the possibility 
that black hole metrics which solve the bulk equations 
with brane boundary conditions, and whose central sin- 
gularities are located on the brane, genuinely correspond 
to quantum corrected (semiclassical) black holes on the 
brane 0, in the spirit of the holographic princi- 
ple [13 and AdS/CFT conjecture 0. 

We recall that our Universe is a codimension one four- 
dimensional hypersurface of vacuum energy density A in 
the BW scenario of Ref. [f| . It is hence useful to introduce 
Gaussian normal coordinates x = (x M ,y), where y is 
the extra-dimensional coordinate such that the brane is 
located at y — (capitol letters run from to 4 and 
Greek letters from to 3). The five-dimensional metric 
can then be expanded near the brane as |17| 

9ab = 9AB\ y=0 + 2 K AB\ y=0 y + £h Kab\ v=0 V 2 

+ -.., (1) 
where Kab is the extrinsic curvature of the brane, and 
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£n the Lie derivative along the unitary four-vector n or- 
thogonal to the brane. We also recall that the junction 
conditions at the brane lead to 

K„u ~ - i (T - A) gitu , (2) 

where T^ v is the stress tensor of the matter localized on 
the brane, and [l7j 

where £ M „ is the projection of the Weyl tensor on the 
brane and f(T)^ v a tensor which depends on T^ v and A. 

The junction conditions in GR [ijj allow (step-like) 
discontinuities in the stress tensor (for example, across 
the surface of a star) keeping the first and second fun- 
damental forms continuous. For thin (Dirac <5-like) sur- 
faces, a step-like discontinuity of the extrinsic curvature 
orthogonal to the surface is also allowed as long as the 
metric remains continuous [19j . Since a brane in RS is 
itself a thin surface, it generates an orthogonal discon- 
tinuity of the extrinsic curvature in five dimensions as 
allowed by GR junction conditions. However, any discon- 
tinuities in the matter stress tensor on the brane would 
induce discontinuities in the extrinsic curvature J5| which 
are tangential to the brane and would therefore appear 
in the five-dimensional metric l[T]l. Such discontinuities 
of the metric are not allowed by the regularity of five- 
dimensional geodesies. Moreover, because of the second 
order term in Eq. Q and considering Eq. Q, we can 
not allow the projected Weyl tensor to be discontinuous 
on the brane either |4l|. One can understand the above 
regularity requirement by considering that, in a micro- 
scopic description of the BW, matter should be smooth 
along the fifth dimension, yet localized on the brane (say, 
within a width of order A -1 / 2 , in units G = c = 1 [2(j )■ In 
any such description, the continuity of five-dimensional 
geodesies must then hold and, in order to build a physical 
model of a star, one has to smooth both the matter stress 
tensor and the projected Weyl tensor across the surface 
of the star along the brane. 

We shall employ the effective four-dimensional (hydro- 
dynamical) equations of Refs. [l7l[T8| in our analysis. In 
general, such equations cannot determine the brane met- 
ric uniquely unless one also knows the bulk geometry. 
However, if the system enjoys enough symmetries, the 
effective four-dimensional equations are closed and can 
thus give some insight about the bulk (i^. I n particu- 
lar, we shall show that, under the simplifying assumption 
that the heat flow is always negligible (no dissipation), 
the knowledge of the full five-dimensional dynamics of 
the most central region of the collapsing star renders 
the whole system "physically closed" when the energy 
density of the star is much smaller than the brane ten- 
sion. By physically closed we mean that the evolution of 
the system is uniquely determined upon further requir- 
ing that the four-dimensional metric become Minkowski 
in the limit of zero energy density and at spatial infin- 
ity (asymptotic flatness). Although considering a non- 
dissipative model appears restrictive, we would like to 



remark that the same kind of models in GR leads to 
paradigmatic examples of black hole formation, beside 
the fact that this is the only case which can be solved 
analytically. 

The Oppenheimer-Snyder (OS) model in GR 
yields the simplest description of how black holes could 
form from collapsing stars. It has however been shown 
that this kind of model is not viable in the RS scenario 
since, if one forces a static exterior outside homogeneous 
and isotropic stars, BW effects produce an effective "en- 
ergy surplus" which is encoded by a positive curvature in 
such an exterior and which cannot be generated by any 
bulk back-reaction. Hence, we expect that this excess en- 
ergy will be released via a mechanism that leads to a loss 
of mass from the star. Although the positive curvature 
in the exterior has no GR description, it can also be ob- 
tained from quantum computations 0] but with the 
opposite sign. The sign mismatch between classical and 
quantum results might be reconciled on recalling that the 
classical anomaly is due to an effective potential energy 
at the boundary of the star which must be released in 
order to have an exterior compatible with the junction 
conditions 23] . It might therefore be possible to change 
the sign of the anomaly just considering that the energy 
surplus should be converted into an effective negative flux 
of energy from the boundary of the star. 

In order to do so, we shall employ a Tolman geome- 
try [24| for the brane star, as it is the only spherically 
symmetric metric which does not allow dissipation of 
energy across the shells of the collapsing star, and the 
Hawking radiation can then be interpreted as the emis- 
sion of gravitons into the bulk. In fact, it turns out that 
the propagation of CFT modes in four dimensions is con- 
sistently described by this mechanism according to the 
AdS/CFT correspondence [l^. Moreover, we require the 
continuity of the Weyl and energy-momentum tensors as 
discussed above, and we shall show that a Tolman brane 
metric corresponds to a general five-dimensional diagonal 
metric with spherically symmetric slices 44] . 

We shall divide the domain of the star into three re- 
gions: I) the "core" , where most of the energy of the star 
is concentrated; II) a "transition region" , which connects 
the core with a tail and finally III) a "tail", where the 
energy density approaches zero. The tail and transition 
region together form an "atmosphere" of the sort that is 
usually employed in numerical simulations of the gravi- 
tational collapse (see, e.g. Ref. [25^. The core is taken 
homogeneous and isotropic (OS-like) for several reasons. 
Firstly, in order to consider a minimal modification with 
respect to the OS model in GR. Secondly, the OS core 
corresponds to an exact five-dimensional solution |29| 
and reproduces the correct Weyl anomaly of quantum 
field theory on the Schwarzschild background |2|, |3| , thus 
making the holographic interpretation clearer. Our main 
results will then be that, in this case, the total energy 
of the system is conserved and that the collapsing 
star "evaporates" until the core experiences a "rebound" 
in the high energy regime (when its energy density is 
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comparable with A) , after which the whole system "anti- 
evaporates" . Moreover, we can find a range of parame- 
ters for which the minimum radius of the collapsing core 
is larger than the AdS length (which sets the scale of 
Quantum Gravity in the BW), thus further supporting 
the qualitative behavior we obtain. 

In Section [HI we shall briefly review a simple holo- 
graphic interpretation of the Hawking radiation in the 
BW. In Section lTTll we shall build a physical model of the 
Tolman type that converges to the OS model in the GR 
limit and show that BW corrections lead to an emission 
of energy from the star and a bounce of the core (see also 
Appendix 0. In Section ITVl we shall analyze in details 
how the (apparent) horizon forms and physical quanti- 
ties related to it, such as the outer trace anomaly, which 
will then be interpreted in terms of the BW corrections 
coming from the brane junction conditions. In section Ivl 
we shall discuss the possible effects of dissipation. We 
shall finally comment on our results in Section I VII 

In the following, we shall use geometrical units with 
G = c = 1 and mostly positive metric signature. 



II. SIMPLE HOLOGRAPHIC PICTURE 

Before trying to "cure" the Weyl anomaly, let us have a 
closer look at the features of the effective energy surplus 
in the exterior of the star discovered in Ref. [a. We first 
note that, as pointed out by different authors |MIT^. l23| . 
the energy surplus reproduces the absolute value of the 
quantum Weyl anomaly computed on a Schwarzschild 
background Q , which is the unique exterior of a spherical 
star in GR. More precisely, in static coordinates one has, 
outside of the star, 



9 M 2 
R» = 

" 2ttA i? 6 



(4) 



where R l \ L \s the Ricci scalar, M the physical mass of 
the star |4(|, and R the Schwarzschild radial coordinate. 
As remarked in Ref. 0], the holographic interpretation 
for such a contribution cannot be of the exact AdS / CFT 
kind, because both classical black holes on the brane and 
semiclassical black holes in GR correspond to strong de- 
viation from AdS and CFT (see also Ref. %). We shall 
indeed show that we cannot reproduce the evaporation 
process if the bulk is simply AdS (i.e. with zero Weyl 
tensor). 

As a first step, we shall show that, if a black hole is 
formed from matter collapsing in the BW, the area of its 
horizon (to first order in A -1 and for a short time after its 
formation) follows the evaporation law for semiclassical 
black holes In particular, we will see that the horizon 
evaporates provided the Weyl contribution is dominant, 
and we may therefore assert that the BW collapse gives, 
to first order in A -1 and for some five-dimensional ge- 
ometries, a good description of the first order quantum 
processes related to it. Let us also note that quantum 



calculations in Refs. 0,0 are performed in adiabatic ap- 
proximation, that is, in some sense, to first order in the 
back-reaction parameter of the quantum theory. 

Following Ref. [||, a unique static geometry which 
matches a collapsing homogeneous and isotropic cloud of 
dust, has a Schwarzschild-like metric with mass function 



(5) 



M = Ms + — m{R) , 
A 



where Mg is the usual ADM contribution (see Section lnTI 
for more details) and 



m(R) = 



3M| 

8 7ri? 3 32irR 



9n 



(6) 



where, in a cosmological background, the constant /i is re- 
lated with the mass of a black hole sitting in the bulk |27| 
and we set the effective four-dimensional cosmological 
constant to zero (since we are just interested in BW ef- 
fects on asymptotically flat branes). 

We denote with R = Ro(t) the radius of the collaps- 
ing object which depends on the proper time r. The 
geodesic equation of motion in the Schwarzschild-like 
space-time determines Rq according to 



2M 2 M s 



Ra 



Ro 4tt\R$ \R$ 



Ml 



(7) 



in which we have selected the case corresponding to zero 
initial velocity for infinite initial radius. We can now 
see how the mass function is changing when the sur- 
face of the star crosses its own horizon, that is at the 
time th when i?o(rn) = Rk = 2Mb = 2A/(th) = 
2 M s + 2m(R B )/X 48] and i? H = -1. Let us define 
the surface area of the evolving "apparent horizon" as 
Aah = 16 7tM 2 (t), for which Eq. JSJ gives 



^ah 
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4A 



Ms 
M 



R 
M 



(8) 



Considering that Mg/M ~ 1 to first order in A 1 , at the 
time t = th, we then have 



4AM, 



H 



or 



9 (i-l 



128 tt A M 2 



(9a) 



(9b) 



The collapse therefore leads to a negative flux of energy 
when the boundary of the star approaches its horizon, as 
expected for the Hawking evaporation, if [/,> 1. 

Since a positive [i generally corresponds to a reinforce- 
ment of the localization of gravity in RS j2^, we can 
assert that an OS region mainly evaporates into gravita- 
tional waves propagating on the brane. We shall however 
see that, for a consistent model of collapsing star with 
continuous density, the sign of the Weyl energy changes 
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from the interior to the exterior of the star, so that the 
evaporation actually ejects energy off the brane via grav- 
itational waves (as suggested in Ref. |T^1)- 

So far, we have not considered any back-reaction on 
the brane metric, and the same flux (|9b|l will reasonably 
be seen by a distant observer for whom d T asymptotically 
becomes a time-like Killing vector. However, the surplus 
energy must be released, since no BW or GR model can 
explain the Weyl anomaly, and this directly implies that 
Eq. (|9bl) probably holds only for a short time about the 
formation of the horizon, as suggested in Q. We shall 
indeed show that this is the case. 



position of the projection of the Weyl tensor on the brane 



1 

8?r 



-Qp. u„ + Q v u^ + , (13) 



corresponding to an effective "dark" radiation on the 
brane with energy density U, pressure hi /3, momentum 
density and anisotropic stress 11^. Note that non- 
local bulk effects can contribute to effective imperfect 
fluid terms even when brane matter is a perfect fluid. 

Bianchi identities supplied by the junction conditions 
produce two kinds of conservation equations 0] : 



III. GRAVITATIONAL COLLAPSE ON THE 
BRANE 

In this section we will study a continuous model for the 
gravitational collapse. In order to see the difference with 
respect to the OS-like model studied in Ref. we con- 
sider a Tolman-like model with a central OS core. The 
star is therefore described as a cloud of dust with falling 
off continuous density and no sharp boundary. The clas- 
sical four-dimensional behavior will be recovered in the 
limit of negligible star density (with respect to the brane 
vacuum energy density A). 



A. General framework 



1. Local conservation equations (LCE): 

p + 6 (p + p) = 

D aP +(p + p) A a = 0: 

2. Non-local conservation equations (NLCE's): 



(14a) 



(14b) 



U + -eU + B a Q a + 2A a Q a +a ab n ab = (15a) 



Q a + 3 © Qa + 3 D Q W + - U A a + D b IU 



+A b IU + a b Q b - Lo b Q b = -E±E 



D Q p ,(15b) 



Following Ref. we can rewrite the BW effective 
four-dimensional Einstein equations with vanishing cos- 
mological constant on the brane as 



G>„ — 8 7r 



Here we have 

->eff cff „. .. i off 7 . off , 



Tiett ett , „ett 1 , ett , eff , i i , 

liV =p u^Uv+p h^v + g ((tl u v) + ir^ , (11) 



(10) 



where D a is the spatially projected derivative (defined by 
D a S b ~... c = //'„//'';...//'', V.. S''' ^, for a = 1, 2, 3), 6 = 
V Q u Q the volume expansion, S a "\..b — u a V Q 5 a "' ...b 
the proper time derivative, A a — u a the acceleration, 
Cab = -D(a u 6) ~ (@/3) h a b the (traceless) shear, and 
u a b = -D [a u b ] the vorticity. 



B. Spherically symmetric dust 



where u M is the unit four- velocity of matter (u^u^ = — 1), 
h^v the space- like metric that projects orthogonally to w M 
{h^u = 9/j.u + Ufj,u v ) and 7r°^ an anisotropic tensor. 

For an isotropic perfect fluid, BW corrections to GR 
are described by the effective quantities 17] 



„cir 



,,cir 



P 1 



_p_ U 

2A + ~p 



7r eff - n 



(12a) 
(12b) 
(12c) 
(12d) 



where p and p are the ( "bare" ) energy density and pres- 
sure of matter. We also employed the following decom- 



For the case with zero pressure {p = 0), that is dust on 
the brane, the quantities in Eqs. (|12afl . I|12b|) and (|12d|) 

reduce to 



P = P 



U 



„cir 



el a 

2A + 3 



Tr eS TT 

v = 11 



(16a) 

(16b) 
(16c) 



Provided the matter density p does not vanish in the 
region of interest, one can use comoving coordinates in 
which u a = (-1,0,0,0). 

In the following, we will only consider the class of five- 
dimensional metrics which are diagonal (sufficiently close 
to the brane at y = 0) and spherically symmetric on the 
brane. In Gaussian normal coordinates, one can always 
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write a bulk metric which is spherically symmetric on the 
brane as 



ds 2 = -N 2 (t, r, y) dr 2 + A 2 (t, r, y) dr 2 

+2 B(t, r, y) dt dr + R 2 (t, r, y) dQ 2 
+dy 2 . 



(17) 



Upon using the restricted freedom to change the four- 
dimensional coordinates on the brane, one can always 
set B(t, r, + ) = [3flj . so that the brane metric reads 



ds 2 



v=o+ 



-N 2 (t, r, 0+ ) dr 2 + A 2 (r, r, 0+ ) dr 2 
-R 2 (T,r,0 + )dn 2 . (18) 



Since we just consider dust as brane matter, from the 
junction conditions at the brane [17j, we also obtain 



= K+(r,r) 



1 dg T1 



2 dy 



v=o+ 



dB 
dy 



= 0. (19) 



y=0+ 



Using the above result together with the bulk symme- 
try Zi with respect to the brane, we have B(t, r, y) = 
y 2 [V(t, r) + 0(y)\. Since the Weyl energy flux is related 
to B by 



Qa 



d 2 B 



dy 2 



(20) 



y=o+ 



one finds that Q a vanishes if V(t, r) — 0, which is in 
fact what we are assuming. The coefficient g Tr then van- 
ishes fast enough on the brane so that, from the five- 
dimensional Einstein equations 



in the limit y 

0= (5) G Tr 



(5) G AB = -Aqab , 

+ , one obtains the condition |4Q 



(21) 



y=o+ 




,(22) 



y=o+ 



where a prime denotes d r and a dot d T . Since our matter 
is pressureless, we can work in the proper time gauge 
N(t, r,0 + ) = 1 and, using the residual gauge freedom 
in defining the radial coordinate r, we obtain 



A(T,r,0+) = R'(T,r,0 + ) 



(23) 



This relation implies a Tolman geometry on the 
brane 



ds 2 = -dr 2 + (R'f dr 2 + R 2 dft 2 



(24) 



where R = R(t, r) is a (generally non-separable) function 
of t and r such that An R 2 (t, r) equals the surface area 
of the shell comoving with dust particles located at the 
coordinate position r at the proper time r. 

With the above symmetries, the vorticity, the accelera- 
tion and the Weyl energy flux vanish, uj a = A a = Q a = 0, 
and we obtain the simplified LCE 



d T p + 6 p = 



(25) 



and NLCE's 



d T u + -eu + a ab u ab = 



- D a U + D b Tl ab 



The volume expansion is also easily computed as 

d T d r (R 3 ) 



= d T [In (R 2 d r R)] = 



d r (R 3 ) 



and for the shear one finds 



1 a h 9 h 

Cab — — T n a b — — li ao 



(26a) 
(26b) 

(27) 
(28) 



where h ao — g ao is the spatial part of the metric l|24|) . 

By symmetry, we expect that the anisotropic pressure 
tensor is diagonal and isotropic in the angular directions. 
Moreover, considering that II Q Q = we have, in such 
adapted coordinates, 



1 



1 



n a 6 = diag - n, -- n, -- n 



and 

& Hah 



■ d T g a b n 



a b 



n 



d T R d T d r R 



R 



d r R 



(29) 



,(30) 



which vanishes in the OS background (homogeneous and 
isotropic space-time) for which 



R(T,r)=g(r)X(T) . 
We then see that the NLCE's become 



U + 



2U 



R!_ 
R' 



U 



n' + 3 %- n 

R 



(31) 



(32a) 



(32b) 



The system of NLCE's is in general not closed, since we 
do not have an evolution equation for LT. However, for a 
sufficiently large physical radius R, the knowledge of LT 
in an extended spatial region together with the asymp- 
totic flatness and the continuity of S^, v make that system 
closed. Let us remark that this also happens in the cos- 
mological perturbative scenario in which one considers 
large-scale evolution of the Weyl tensor plj . 

The LCE (|25|l integrated over the spatial volume 
Vh dr d9 d(f> = sin 0/3 d r (R 3 ) dr d0 d</> implies 

d T m p = , (33) 

where we have introduced the "bare" mass function 

m p (r) = / p(r,x)d x (R 3 (t,x)) dx , (34) 
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or, equivalently, 



p{r,r) 



AttR 2 R' 



(35) 



The meaning of Eq. I|33|) is that, since we have chosen 
a comoving reference frame and p = 0, the "bare" en- 
ergy contained within a sphere of fixed coordinate radius 
r cannot change in time, although the physical radius 
R(t, r) of such a sphere decreases during the collapse. 
We can now consider the (r, r) Einstein equation, 

GT ^ = -w = - 8 ^ off ' (36) 

which yields the equation of motion 

k , 2M(^) + f ( T) i 

R{T,r) 

where we have introduced the "effective" mass 



M (t, r) 



4"7T 

~3~ 



p cB {T,x)d x (R 3 (T,x)) dx . (38) 



Since we want a flat brane for M = [moreover, the 
center of the star is at rest, R(t, 0) = 0], it must be 
F(t) — and we finally obtain 



R 2 (r,r) 



2M(t, r) 
R(r,r) 



(39) 



Let us note that the effective mass is not constant in 
general. In fact, 



M(r,r) = ^-8 T J 



cir 



d x (R s )dx 

f x+ u)a x (R>) 



dx . (40) 



For the particular case II = 0, one then obtains 

M(r,r) = jf + |) d x d T (R 3 ) dx , (41) 

where we have used both the LCE and the first NLCE. 

A very important result which follows from the LCE 
and NLCE's is that, if the brane metric is asymptotically 
flat, the anisotropic stress LT ^ whenever p' ^ 0. We 
can prove it by showing that LT = is not compatible 
with asymptotic flatness and the LCE and NLCE's. On 
combining Eq. (|32a[) with Eq. 125|) for LT = 0, we obtain 



U = U(r) p 4/3 



(42) 



where U(r) is a time-independent integration function. 
From Eq. I|32b() with LT = 0, one instead obtains 



U 



2A 



+ F(r) , 



(43) 



Core Transition region Tail 

FIG. 1: Density profile. 



F(t) being a spatially-constant integration function. 
Asymptotic flatness requires that 



lim U(r,r) — lim p(r,r) = , Vr 



(44) 



which implies F(t) — 0. On now combining Eq. I|42|) 
with Eq. (|43|l . we get the relation 



U(r) 



3p 



2/3 



2A 



(45) 



which obviously contradicts the assumption p' ^ 0. This 
implies that, for a continuous distribution of dust for 
which p' ^ 0, one must have II ^ 0. This result supports 
the holographic interpretation as it resembles very much 
a property of the renormalized quantum stress tensor on 
the Schwarzschild background @ . 



C. The model 

As mentioned before, we shall divide the star in three 
regions [HJl (see Fig. for a qualitative picture): 

I : the "core" (0 < r < ro), where p' = and one has 

an OS [2^ behavior; 

II : the "transition region" (ro < r < r s ), with Tolman 

|24| behavior due to p p' being non-negligible; 

III : the "tail" (r > r s ), with p p' ~ 0. 
Moreover, we define the dimcnsionlcss parameter 



Po/A , 



(46) 



where po = p( T — 0, r = 0) is the initial core density. 
Such a parameter is assumed small, since the system 
is initially in a low energy regime (from the BW point 
of view) and relevant quantities can thus be expanded 
in powers of e for sufficiently short times (or sufficiently 
large distance from the core). 

A basic feature of both Tolman and OS models in four- 
dimensional GR is that the bare mass function at fixed 
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comoving radius is constant in time and remains well 
defined during all the collapse. Therefore, dust shells 
of different comoving radius move along geodesies solely 
determined by the inner geometry and reach the central 
singularity (R = 0) at increasing proper times (Tolman 
model) or at the same proper time (OS model). In the 
former case one can have an enlarging apparent hori- 
zon |52j, while in the latter just an event horizon forms 
at the star surface |3(J. 

In the BW, the role of the bare mass is taken by the 
effective mass M of Eq. (|38|l . which will be shown to 
diverge whenever R — * 0, thus making the whole four- 
dimensional space-time singular. To avoid this case, 
which is mathematically admissible but physically un- 
likely, one has to include a sufficiently negative contribu- 
tion to the mass coming from the projected Weyl tensor. 
As we discussed in Section^ this will generate an Hawk- 
ing flux near the forming horizon, and we shall further 
show that the effective mass completely evaporates at a 
finite star radius, after which the collapse changes to a 
re-expansion (or "anti-evaporation" process). This case 
of BW collapse and rebound cannot be related to the GR 
behavior perturbatively (in e ~ A -1 ), since none of the 
shells reach R = 0, but we incidentally note that it seems 
in agreement with the uncertainty principle of quantum 
mechanics jH^]. In fact, a "bounce" in the trajectories of 
the collapsing matter caused by quantum gravitational 
fluctuations had already been found in an improved semi- 
classical analysis of the OS model [32| . 



1. The core 

We first recall that the bulk solution which corresponds 
to the OS core of the star is perfectly regular in five di- 
mensions far from the space-time singularity p9l |. Fur- 
ther, since p' = 0, the system of relevant equations is now 
closed. In fact, we have II = 0, and the NLCE's reduce 
to the one equation 



effective mass l|38|) is then given inside the core by 



U+ -QU = , 



which is solved by 



U = - 



27pr 4 e 



128 7T 2 po i? 4 



(47) 



(48) 



where p is a constant. 

The physical radius R can in general be written in 
the factorized form (|31|l and the coordinate r can be so 
chosen that 



9 \ 1/3 r 



(49) 



M(r,r) = M s 



+ 



r 
ro 
9 c 



32 it po 



(2 M s ) 
3i? 3 



E 

R 



(50) 



where the first term in the r.h.s. is the usual bare mass 
and the remainder represents the BW correction. 

The above effective mass would diverge for R — > (this 
also occurs for a general Tolman core, see Appendix |A"|) . 
The point R — is the usual central singularity, which is 
harmless (at least when covered by an horizon) in four- 
dimensional GR, since the bare mass is constant and fi- 
nite. In the present case, however, the diverging effective 
energy makes the whole space-time singular. In order to 
see this, let T be the proper time at which the OS core 
hits the singularity. From the equation of motion (I39"|) 
one has 



RR Z 



r>r a 



= 2M(T,r ) +8tt / p cS R 2 R' dx . (51) 



in which Ms is the total bare mass of the OS core. The 



Since M(t — * T, ro) — > oo, either the second term in the 
r.h.s. is finite and the total effective mass diverges at any 
r > ro, thus making the whole exterior singular, or it 
equals — M(T, ro) + f(r), with f(r) a regular function, 
in order to compensate for the diverging core energy. In 
the latter case, the Weyl energy becomes everywhere in- 
finitely large and negative and, since G a a ~ U, the whole 
four-dimensional Einstein tensor is singular. Although 
such singular evolutions appear mathematically allowed 
by the equations, in the following we shall not consider 
them since, from the BW point of view, either they pre- 
dict a catastrophic end of the Universe induced by astro- 
physical events or, more reasonably, they suggest that a 
more fundamental model must be used. However, in the 
latter case we expect for large black holes that a huge 
energy flux would be emitted towards infinity well before 
the OS boundary approaches the Planck length. This, of 
course, would be ruled out by astronomical observations. 
From the holographic perspective we are interested in 
here, only the non-singular solution is relevant. In fact, 
it is only in this case that the star continuously "evapo- 
rates" (before the bouncing) by emitting a Hawking-likc 
energy flux at the moment when the OS horizon forms, as 
we shall see later. Moreover the bouncing solution seems 
to be compatible with some proposal for the quantum 
black hole formation 

In order to avoid the singular cases, one must have p 
positive and large enough so that each shell will bounce 
back after reaching a minimum radius where the corre- 
sponding effective mass vanishes 54] . The Weyl tensor, 
holographically interpreted as the quantum back-reaction 
on the brane metric (see j2^| and References therein), 
then contributes the "evaporation" term proportional to 
p in Eq. (|50|l , which dominates at relatively low energies; 
whereas, the BW correction to the matter stress tensor, 
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X 



FIG. 2: Qualitative behavior of the shells and core potential 
for p > p c (solid line) and p < p c (dashed line). For p = p c 
(dotted line), the peak of V equals the shells energy E — 0. 



holographically interpreted as the matter quantum trace 
anomaly 01 > yields the "anti-evaporating" term propor- 
tional to Mg in Eq. (|50|) , which increases with the energy. 

Upon inserting the effective mass (|50|) in the equation 
of motion (|39[) . one obtains an equation for X(t), 



X 2 = 



6 l/3 



cp 



9X 277rpo^ 4 
-V(X) , 



24 ir po Mg /3 X 2 



(52) 



in which there is no dependence on r. This shows that 
the system remains "rigid" through the bounce: no shell 
crossing occurs and all shells reach their minimum radius 
at the same proper time. Like for the classical OS model, 
it is thus sufficient to consider the evolution of the core 
surface at r = r$ and we correspondingly define Mq(t) = 
M(r,r ) and R {t) = R(r,r ). 

From the form of the potential V in Eq. I|52|l (see also 
Fig. [21, one can see that the term proportional to p 
behaves as a repulsive (angular momentum-like) force 
and the bounce occurs whenever there is a positive peak 
(since the energy of collapsing shells E = for our choice 
of initial conditions). There will in general exist a critical 
value p c — n c (Ms, Po, e) such that one has the bounce 
for p > p c , otherwise Rq — > and Mq diverges. For 
p = p c the two turning points of the potential coincide 
and the shells would take an infinite proper time to reach 
the minimum radius (of course, this would only occur if 
one neglected any perturbations, and we shall not fur- 
ther consider this special case). In Fig. |3|we display a 
typical trajectory of Rq, along with the corresponding 
time evolution of the core effective mass Mq, for p < p c 
in panel (a) and for p > p c in (b). In the latter case, 
after the core surface has reached the point of zero ef- 
fective mass, it will bounce back transforming the whole 
collapse into an "explosion" , which evolves as in (b) with 
the time reversed. Although we are not able to describe 
the dynamics of the atmosphere at very high energies 
(e.g. around the bounce) by means of our perturbative 
analysis, on considering that the continuity of the total 
energy-momentum tensor would be spoiled if the shells 
crossed (55j , one finds that all the shells (both in the core 




00 

FIG. 3: Typical evolutions of the core radius Rq(t)/Ro(Q) 
(solid line) and effective mass Mq(t)/Mo(0) (dashed line) for 
p < p c (a) and for p > p c (b). Units are arbitrary. 



and the atmosphere) must begin to re-expand. Because 
of energy conservation at infinity, this "reversal of mo- 
tion" in the atmosphere would generate an instantaneous 
distributional (Dirac <5-like) term in the Ricci scalar, as 
it was also found in a semiclassical treatment of bounc- 
ing solutions [40| . Such a singularity in turn means that 
a detailed description of the collisions between matter 
shells at fixed r inside the atmosphere and with the core 
must be taken into account at that point. A microscopic 
description of the shells goes however beyond the scope 
of the present paper and we just wish to make a remark. 
In practice, during the bounce the collisionless descrip- 
tion of dust must be relaxed by introducing an effec- 
tive short distance potential which results in an effective 
equation of state for the atmosphere. Since our system 
is non-dissipative by construction, the scatterings should 
be completely elastic and the equation of state of the 
polytropic type (see, e.g. Ref. (25[). Of course, this prop- 
erty can be viewed as an artifact of our simplified model, 
whereas in a more realistic situation some energy will be 
dissipated from both the core and the atmosphere, as we 
shall discuss in Section Ivl 

The turning points of Rq can be found analytically by 
solving the cubic equation V = 0, but their expression 
is rather cumbersome. It is instead easy to determine p c 
exactly by observing that the peak V — is located at 
X = X p = (4 Af s /3) 2/3 p~ 1/2 for some values of M s and 
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fx, and, in general, 



V(X p ) = y/ji- 



3 e /1 3 / 2 - 32 7T po M 2 



48 6 1 /3 7rpo M s s 



8/3 



(53) 



Hence, V has two positive zeros if and only if V(X P ) > 0, 
that is when 



li > 



{ 32irp a 
V 3e 



Ml 



2/3 



(54) 



We remark that the bouncing is a high energy effect 
compared to e (since fi c ~ e~ 2 / 3 ), whereas the evapora- 
tion also occurs at low energy. In fact, the core effective 
mass is given by 



M (t) = M s + 



9 c 



32 ir po 
and its time derivative is 



(2M S ) ; 

3i? 3 



M (t)=- t 



9e 



32 7r po 



2 M s 



i?0 



(55) 



1* • (56) 



/? 2 



Recalling that Rq < during the collapse and that 
i?o(0) ^> 2 Ms (the initial core radius must be outside 
the GR horizon), we see that the evaporation sets out 
at the beginning of the collapse when the star is still in 
the low energy regime. Moreover, thanks to the condi- 
tion (|54p . one can easily show that Mo < 0, at least until 
the radius bounces back. 

In particular, the minimal radius R m in is given by 



> x- 1 ' 2 



M s 
A-V2 



1/3 



(57) 



The holographic description is expected to hold only if 
the AdS length I ~ A -1 / 2 is much shorter then the typ- 
ical lengths of the process we are considering 15]. The 
shortest length in our system is obviously given by i? m i n , 
for which we should therefore have 



where the minus sign (Mo < 0) holds during the collapse 
and the plus sign (M > 0) after the bounce, and Ro(t) 
can be determined to zeroth order in e. 

From now on, we shall just analyze the collapse, since 
the explosion is the time reversal of the latter in our case 
where there is no dissipation. Since it is the core which 
first enters a high energy regime, we can obtain a (rather 
conservative) estimate for the error made if we truncate 
expressions to first order in e by comparing Mo to first 
and second order by means of the function |56( 



A(r) 



dlMo 



ft Mr 



e=0 



e=0 



(61) 



and consider that our approximation is good if A < 0.5. 
The analytic expression of A is extremely involved and 
we just show a few plots in Appendix from which the 
dependence on Ms and /i can be qualitatively inferred. 

A physical upper bound on \p\ can be placed by con- 
sidering that the BW correction to the core bare mass 
for astrophysical objects must be much smaller then the 
bare mass, 



|M -Ms| «M S 



(62) 



(at least) until the core approaches the GR horizon (Ro 
2 M s ), and Eq. then yields 



i i „ 6471-po ,, 2 
W < 9e M s 



(63) 



For astrophysical objects one also expects AMI 3> 1, so 
that p, a ^> fi c . Moreover, the limit (|63|l assures that the 
formation of the OS (apparent) horizon, occurs before 
the bouncing. However, this upper bound cannot likely 
be used for small black holes for which we expect a strong 
Hawking evaporation even at the formation of the first 
horizon. 



2. The transition region 



Rn 



» 1 



A-V2 

From Eq. I|57|l . we thus need 

M s > A" 1 / 2 



(58) 



(59) 



that is, the Schwarzschild radius of the star must be much 
larger than the AdS length as one would have expected. 
Furthermore, from Eq. (|59|l it also follows that p c ^S> 1. 

In light of this remark, in the following we will study 
the evolution of the whole system to first order in e, to 
which we have 



M (r) = T 



9 c 



32 ir po 



4M| 



Rlir) 



V2Ms~ 
R 5 /2 (r) 



(60) 



For ro < r < r s , we are in the transition between 
two regions of almost constant density. Since in the GR 
model p = for r > r , the energy outside the OS star 
is entirely a BW correction. The density therefore must 
decrease rapidly from a value which is of order e° to a 
value of order e. This can be formalized as 



4 it f r 

m p (r;r ) = — / p (R 3 )' dx = O 



(e) , (64) 



where ro is again the border between the regions I and II 
and the LCE as usual guarantees that m p remains con- 
stant. Moreover, since the transition is overall a BW 
effect, we can take 



r s - r = OiX- 1 ' 2 ) = 0(e) 



(65) 
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and therefore 



R( T ,r)-Ra(T)=0(e) 
Since U = 0(e), we also have that 



(66) 



m u (T,r;r ) 



47T 



U (R 3 )' dx = O (e 2 ) , (67) 



for ro < r < r s , and the contribution of U to the effective 
mass in region II can be neglected. Although in the tran- 
sition region we have no control on the projected Weyl 
tensor, we can still regard the system as closed since the 
Weyl contribution does not affect the evolution at the 
level of approximation we are considering. Combining 
these results, we obtain that, to first order in e, the ef- 
fective mass is given by 

M(t, r) ~ Mo(r) + m p (r; ro) + m u (r, r; ro) 

~ M (t) +m p (r;ro) . (68) 

This implies that, to first order in e, 

M(r,r) ~ M {t) . (69) 

for r < r < r s , in agreement with the condition (|65|l . 
and we can conclude that, since M < at r — ro, it will 
remain negative (and substantially unaffected) through- 
out the border of the transition region r = r s . 

3. The tail 



As in the transition region, p = 0(e) for r s < r, and 
TO p( r ! r s) = 0(e). Furthermore, we can now consider that 
in this regime p' p/X = O (e 2 ), so that bulk gravitons are 
decoupled from brane matter. The Weyl contribution is 
however of the same order, 

m u (r,r;r s ) = ^ J* U (i? 3 )' dx = 0(e) . (70) 

We recall that the effective Einstein equations imply 



that the Ricci scalar R^,, 



!7rT eff , that is 



R^ = ^d r [Rd^]=^[p- P - 



(71) 



Upon integrating over regions II and III and taking into 
account Eq. (|64[) . we thus obtain, to first order in e, 



RdlR 2 



2m p {r;r ) , 



(72) 



for ro < r. From the equation of motion l|39|l . the above 
relation yields 



5m 

R 



RdlR 2 



2m p (r;r s ) - 2to^(t, r;r s ) 



-2m u (T,r;r s ) 



(73) 



now for r s < r. On further considering Eq. (|67|l . we 
obtain 



R ■ 
— M 
R 



R ■ 
— M 
R 



-2mu(r,r;r s ) 



(74) 



Since Mq = 0(e), we can use the zeroth order equation 
of motion for the shells at fixed r > r s , 



R 2 (r,r) 



2 M s 
R(r,r) ' 



(75) 



where Ms is again the total bare mass of the OS core. 
Solutions to the above equations can be written as 



R(r,r) 



1/3 



[/(r)+T-r] 2/3 , (76) 



where the function f(r) is monotonically increasing in r 
and such that R(t, r) is continuous across r = r s . There 
is no loss of generality in assuming that f(r) — r — c with 
c a constant, since changing / is tantamount to redefining 
the coordinate r. In particular, on considering Eq. 16611 . 
we can set c = ro to zeroth order in e (for a discussion of 
T, see Appendix IbI). 

One can now prove a general result which holds irre- 
spective of the specific solutions for U and II. Since, for 
r > r s , 

M(t, r) ~ M (t) + m p {r; r ) + m w (r, r; r s ) , (77) 

one has that 

M(t, r) ~ M (r) + m u (r, r; r s ) , (78) 

and, for any given R = R(t, r), Eq. (|74|l becomes a differ- 
ential equation for mu — mu(T, r; r s ), whose form further 
simplifies on taking into account the approximation 17511 , 



2 M s 

mu + 2\j -^g- m u = M 



Rq 
R 



3/2 



- 1 



(79) 



Since R(r,r) > Rq{t) for r > r , Eq. implies that 
mu cannot remain zero in the tail (note that the r.h.s. is 
positive for M < 0). 

In order to proceed, we now assume that: 

(i) linv^co R(r,r) = oo |H3' an d 

(ii) the effective mass be finite at spatial infinity, 

lim M(t, r) < oo , Vr > . (80) 

r — >oo 

Since Mo(r) always remains finite if p. > p, c (i.e. when 
there is a bounce), and the bare mass of the tail is finite 
(and small) by construction, this implies that 



lim mu (r, r ; r s ) < oo , V r > 



(81) 



Since asymptotic flatness ensures that at large distance 
from the core the low energy approximation holds, we 
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FIG. 4: Typical contribution of U to the effective mass at five 
subsequent times (from horizontal axis for r = to upper 
curve for r = T/20) for e = 10" 4 , M s = p = 1, T = 10 and 
Ai = 5000 > 4824 = fi c . 



R/M s 




M s = po = 1, T = 10 and /z = 5000 > 4824 = /i c , 
but we note that different values of these parameters do 
not change the qualitative behavior of mu and M. In 
the following, we shall use these values of the parameters 
for all the numerical computations and plots. They can 
in fact be considered as the case of a small black hole 
for which, however, the Holographic bound is satisfied. 
The only purpose of the plots is to show more clearly the 
qualitative behavior of the processes involved as well as 
to support our perturbative expansion for any physical 
interesting cases. 

Since mu(r,r;r s ) increases monotonically in r (start- 
ing from zero at r — r s ) and, because of Eq. JHHJl, the 
effective mass of the star outside the core will decrease 
by releasing gravitational waves off the brane and into the 
bulk |28( . After the bounce, since the core will eventually 
re-enter a low energy regime in which Mq is the same as 
the one found before with opposite time evolution, we ex- 
pect that mu will also evolve backwards so as to ensure 
the general condition l|83|) . At a time equal to twice the 
time of the bounce, we should therefore have mu = 0, 
corresponding to the initial state with zero Weyl energy. 
This behavior is related to the non-dissipative nature of 
our model, and we shall later discuss the possible effects 
of dissipation. 

IV. BLACK HOLE FORMATION AND 
EVAPORATION 



FIG. 5: Time derivative of the effective mass M(t, r) ver- 
sus the physical radius R(t, r) for four values of r = 
T/100, . . . , T/20 (upper curve to lower curve) and the same 
parameters used in Fig. 2] 

can take the limit r — > oo (equivalent to R — > oo at fixed 
time) in Eq. (|Tf)|) and finally obtain 

lim m u {T, r; r s ) = -Mq(t) , Vr > , (82) 

r — >oo 

or, from Eq. (JJSJ, 

lim M(r,r) = , Vr > . (83) 

r — >oc 

To summarize, we have shown that if the total effective 
mass at spatial infinity is finite at the initial time r = 
0, it will always remain constant (for a bouncing core 
evolution with /! > /x c ), so that the total effective mass 
of the collapsing dust star is actually conserved. 

Eq. i|79|l can be solved analytically with a generic ini- 
tial condition m^(0, r; r s ). It is particularly interesting 
to consider the case mu{0,r;r s ) — so that there is ini- 
tially no energy stored in the Weyl component U (the 
bulk metric is AdS at low energies). In this case, us- 
ing Eq. f75)l. we have that mu > during the collapse. 
This yields the curves of Fig. 0] and the time derivative 
of the effective mass of Fig. El m which we set e = 10~ 4 , 



We now proceed to analyze the model developed in the 
previous Section near the (forming) horizon. 

A. Horizons 

We recall that shells of constant r reach the (apparent) 
horizon at the time r = th (r) when 

Rfa,r)=-l, (84) 

provided R(tb,x) > —1 for a; > r (at least locally). 
From the equation of motion i|39fl . this is equivalent to 
R(rn,r) = 2M(rn,r), which means that surfaces of dif- 
ferent comoving coordinate r reach the null surface (hori- 
zon) at different proper times. We can equivalently define 
r = t*h (t) as the value of r at which the horizon is formed 
at the time r. 

The function th(t) is of course model dependent and 
affects how the effective mass evaluated on the horizon, 
Mh(t) = M(t, rn(r)), changes in time. In fact, its total 
time derivative contains two contributions, 

^5=M H + M^r H . (85) 
clr 

The first term in the r.h.s. originates from the intrinsic 
time dependence of the BW effective mass at constant r 
which we have studied in the previous Section, is a first 
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order effect in e and would vanish in GR. The second 
term accounts for the mass change due to the (possibly) 
variable number of shells included within the horizon and 
depends on the detailed form of the atmosphere. Since 
we have assumed that our model is OS to zeroth order in 
e, we have Mjj = 0(e) outside the modified OS boundary 
(r > r s ~ r ). 



d&f -0.02- 

"37 



In the core 



Since in this region the model is OS, the velocity 
|i?(r, r)| increases monotonically in r at fixed r. There 
is therefore only an (apparent) horizon at the boundary 
r = r when Rq(t) satisfies Eq. I|84|) 58|. This occurs at 
the time (to zeroth order in e) 



t 



FIG. 6: Total time derivative of the effective mass Mh(t) 
versus the time t = (t - t§ s )/T for the BW model (solid 
line) and for the Hawking law (dotted line) for e — 10~ 4 , 
Ms = po = 1, T = 10 and p, = 5000 > 4824 = p c . 



t h )S = T--M s 



(86) 



where T fixes the time scale of the collapse (this would 
be the time at which the star hits the central singularity 
in GR; see also Appendix [5J . On the event horizon of 
region I, we then get 



dM H 
dr 



Mu 



Mo(t£ & ) 



9 (M-l) e 
128 7tpo-Mh 



(87) 



which is precisely the Hawking flux obtained in SectionlTTI 
once we replace the definition e = po A -1 . 



2. In the transition region 

Inside this Tolman region, the horizon for a given shell 
will be reached at proper time r > t§ . The partial time 
derivative of the effective mass on the horizon will then 
scale according to Eq. I|69|) . In particular, on considering 
again that r — r — 0(e), the total derivative scales as 



dA/ H 
dr 



Mk{t,tk{t)) ~ M (r) ~ M (t§ s ) , 



in which the last approximate equality follows from the 
condition (|66[1 . This implies that the flux at the OS hori- 
zon will continue up to the time when ru(r) = r s . 



3. In the tail 

This is the most interesting part, since the above re- 
sults for the transition region allow us to approximate the 
boundary of the modified OS star as the sphere r = r s . 

Since both Mh and in Eq. (|85|l are of order e, we 
can use the zeroth order Eq. I|75[) in order to determine 
the evolution of the horizon which therefore stays at 



R(T,r H (T)) ~ 2 M s 



(89) 



or r H (r) 



^os 



Eq. H85|) with fn — 1, we then find 
dM H 



Upon inserting Eq. I|77|) into 



dr 



M (r) + m w (r, r H ;r s ) 
+m' p (r H ; r s ) + m' u (r, r H ; r s ) , (90) 



for r > t§ s 



As opposed to the other terms in Eq. I|9U|I . 
the contribution given by M^ depends on the specific 
profile chosen for p, and is present in GR Tolman model 
as well. The increase of the mass at the horizon induced 
by this term is simply due to a flux of matter flowing 
towards the center of the star which makes the appar- 
ent horizon grow. Obviously m' p is positive, does not 
explicitly depend on time [but just via rn = J"h(t)] and 
decreases for increasing ru (or, equivalently, for increas- 
ing time t). The deviation of the smooth energy density 
of the atmosphere from the OS outer vacuum should be 
local, hence very much concentrated near the OS bound- 
ary. This implies that the profile of the density should 
decay very fast 59] . With this in mind, one can choose r s 



in such a way that p ~ O(e ) and r 



O(e), so that 



m' p is negligible to first order in e. We will then not con- 
sider its contribution to the total variation of the mass at 
this stage. The term m! u is instead determined uniquely 
by Eq. I|79|l and the initial condition for U (which we 
naturally took as zero Weyl energy). 

From Eqs. I|79|l and l|89|) . the first contribution is easily 
approximated as 



M H ~ M (r) 



Ro(t)\ 3/2 m w (r,r H ;r s ) 



2 Ms 



M s 



9 c 



64 7T po M s 



4MJ 

OS. 



Ro(r) 



mu(T, r +T - r H )b ;r s ) 
M s 



(91) 



in which we finally used Eq. I|6l)|l and i?o( T ) is of the form 
l |7B) with f(r) = 0. 
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FIG. 7: Error for the solid line in Fig.© Note that t < 
correspond to times before the horizon formation (r < Th S ). 
The dashed line marks the limit of validity A < 0.5. 



Since the analytic expression for rnu is too compli- 
cated to display, we compute the total time derivative of 
the effective mass at the horizon from the solutions of 
Eq. H79fl numerically. For the values of the parameters 
used in Figs. 01 and [5] the result is plotted in Fig. [fj] up 
to the time when the error estimate A ~ 0.5 (see Fig- EJ - 
We can see that the flux is smaller with respect to that 
predicted by Hawking, and this behavior remains for dif- 
ferent values of the parameters. In particular, decreasing 
e or /i, as well as increasing Mg, reduces the luminos- 
ity, as expected, and keeps our approximation reliable 
for longer times. 

The conclusion is that, although the evaporation sets 
out according to Hawking's law, the back-reaction on the 
brane metric subsequently reduces the emission until the 
effective mass of the core vanishes and its radius bounces 
back with Mq that becomes positive |60|, There will be 
an interval of time during which the error A is large and 
our first order analysis outside of the core breaks down. 
However, after a finite amount of proper time, A will 
become small again and the system will evolve back to the 
initial condition through a sequence of states obtained by 
inverting the time in the above solution. 



B. Luminosity 

A distant observer experiences an impinging flux of 
energy during the collapse, whose total amount must be 
calculated from the horizon to infinity (since the region 
inside the horizon is causally disconnected from such an 
observer). Since the total energy from the origin to in- 
finity is conserved and the process extracts energy from 
the hole, we expect that the measured flux is positive. 

After the horizon has formed on the boundary of the 
star (more explicitly, for 7"h > T s ~ tq and r > i~§ s ), one 



has 



** 55 d7 
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lim — 



„cir 



R 6 ) dr 



)"h(t) 



= li m M(,0-^ = -^,(92) 
dr 



dr 



in which we used the conservation of the total effective 
mass (|83t . Further, since d T — > dt for r — > oo one finally 
obtains the luminosity 



dM H 
dr 



(93) 



The flux seen by a distant observer therefore shows the 
same dependence on the mass Mr as the semiclassical 
expression when the horizon is first forming, and subse- 
quently decreases to zero (before it becomes negative). 
However, since this happens after the apparent horizon 
begins to form, a distant observer might have to wait an 
infinite amount of time to measure a vanishing flux. 

We now consider the only instant when the Hawking 
radiation actually equals the BW result, that is at the 
OS boundary when rn(r) = r*o. Reintroducing the New- 
ton's constant G in units c = H = 1 (so that M has the 
dimension of a mass in this Section, and not of a length) 
and the definition of e, we obtain 



9fc-l) 
128 ttG 4 AAf2 



(94) 



which we can compare with the semiclassical luminosity 
as calculated in the Schwarzschild background |3J| 



a 



G 2 M 2 { 



(95) 



where a is a dimensionless coefficient which depend on 
the quantum field theory chosen. 

An astrophysical object has M ^> 10 31 GeV and we 
can therefore use the result (|9"5*|) of Ref. , that is a = 
q N 2 ~ 7.74- 1CT 3 A^ 2 , where N is the number of particle 
species appearing in the quantum theory 16 1| . On now 

using the lower (/i c ) and upper (/i a ~ nV 2 ) bounds for 
/x as given in Eqs. I|54|l and (|62fl . we obtain a limit on 
the number of species that can take part in the Hawking 
process, 



128 7ra V 3 



32 



2/3 



M 4 



1/3 



< N 2 < 



G M 2 
2a 



(96) 



The result of Ref. [33| must be valid for any mass M ^> 
10 31 GeV Therefore, for the upper bound we can safely 
consider M ~ 10 35 GeV, so that 



10 43 ( 



GeV 4 
A 



1/3 



< JV < 10 



X', 



(97) 



For consistency, we must also have 



A > 10 30 GeV 
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Considering that the AdS length i = G A must be 

much larger than the Planck length, we finally obtain 

~ 32 mm < £ < 10" 9 mm , 



10" 



(99) 



This bound for the AdS length is three orders of mag- 
nitude better than the best constraint found in Ref. PJ| 
considering the time scale of primordial black hole evap- 
orations. 



C. Trace anomaly 

Strictly speaking, there is no trace anomaly in our ap- 
proach, since we have included the back-reaction of the 
effective matter on the brane metric. However, in order 
to compare with known results without the back-reaction, 
we can define the trace anomaly 1Z as the sum of the Ricci 
scalar and the trace of the bare stress tensor 62] . From 
the effective Einstein equation i|71[l one readily obtains 
(see also Ref. |3^) 



1Z 



1 



,4 



2vrA \R 2 R 



At the OS boundary, r = rp, we then have 



H 



9 Ml 



(100) 



(101) 



2ttA R 6 ' 

which is the quantum Ricci anomaly of Ref. Q with 
the correct sign at the collapsing boundary. It is then 
clear that the sign mismatch found in Ref. Q was due to 
the choice of a non-smooth energy density and that, by 
adding a tail, we have described how the excess energy 
stored in the OS boundary is released. 

For r > ?*o, we have 1Z = 0(e 2 ), which is therefore 
negligible from the point of view of our analysis be- 
cause, as we have shown, the Hawking flux diverges in 
time whereas the BW one remains finite outside the OS 
boundary. How the back-reaction on the brane metric 
gradually annihilates the Ricci anomaly in the transition 
region can only be understood by introducing specific 
models for the tail which must also be consistent with 
the five-dimensional problem, and goes beyond the inves- 
tigation we want to present here. In any case, if the holo- 
graphic analogy holds, the anomaly of Ref. [2j must just 
be effective at the boundary of the OS core and decrease 
to zero at the modified boundary of the star (r = r s ). 
It is so because the back-reaction on the brane metric 
must be consistent with the modified Einstein equations, 
whereas in Ref. the Einstein equations are just solved 
for the background. 

V. DISSIPATIVE COLLAPSE 

We now wish to discuss the possible effects of a dissi- 
pative term of the form Q a =/= in our model, although 



including such an energy flow would require heavy nu- 
merical investigations of the full five-dimensional equa- 
tions and goes beyond the scope of the present paper. 

Since, the OS core of our dust star is non-dissipative by 
construction, the off-diagonal term V(t, r) — for r < tq 
in the bulk metric l|17[l . as we explained in Section lIlIBI 
This implies that one can have Q a (T,r) ~ V(r,r) ^ 
only for r > tq with an OS-like core. In this case, com- 
patibility with the Hawking effect would constrain the 
heat to flow from the OS boundary towards infinity, and 
the energy of the transition region and tail would thus 
be dissipated away completely after a suitable amount of 
time. In the meanwhile, the core should keep bouncing 
back and forth between its initial condition and the state 
with vanishing effective mass, since its evolution cannot 
be affected by Q a in the atmosphere. The net final result 
should thus be that the system converges to the model 
with an empty exterior discussed in Ref. |6J , which we al- 
ready know is not acceptable in the BW. This argument 
shows that a non-dissipative core is most likely incom- 
patible with a heat flow in the external region and that 
the condition Q a = should therefore not represent a 
real restriction for the model we have analyzed. 

Of course, a more realistic model for a collapsing star 
should also have a dissipative core and one should con- 
sider a non-vanishing Q a everywhere, as well as a non- 
vanishing flow of matter. The (absolute value of the) to- 
tal (holographic) flux of energy measured far away from 
the core would then be larger than the one from a non- 
dissipative core, and therefore closer in value to the four- 
dimensional Hawking flux. We nonetheless expect that 
a global horizon does not form, since the total outgoing 
flow will make the star "evaporate" until all the initial 
energy has been radiated away. This can happen either 
before or after the bouncing, which does no more allow 
the star to come back to its initial condition because of 
dissipation. In fact, we expect that no singularity forms 
even in the general case because the badly diverging part 
of the Ricci scalar proportional to the squared energy 
density in Eq. (|100(l . which arises from the junction con- 
ditions on the brane, cannot be canceled by the Weyl con- 
tribution. The singularity must then be avoided either 
with the help of the Weyl tensor or by means of severe 
modifications to the matter profile due to BW effects. In 
the former case, the star will still bounce, whereas in the 
latter it will completely "evaporate" before reaching the 
singularity. This anyways remains an open question that 
cannot be addressed here. 



VI. CONCLUSIONS 

Inspired by the conjecture that classical black holes 
in the BW may reproduce the semiclassical behavior of 
four-dimensional black holes, we have studied the grav- 
itational collapse of a spherical star of dust in the RS 
scenario in order to clarify the underlying dynamics that 
leads to this interpretation. Regularity of the bulk geom- 
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etry requires continuity of the matter stress tensor on the 
brane and can lead to a loss of mass from the boundary 
of the star. We have in particular shown that, excluding 
energy fluxes coming from the bulk Weyl tensor, a col- 
lapsing spherical star must have a spatially anisotropic, 
although isotropic in the angular directions, atmosphere, 
in order to have asymptotically flat solutions. Interest- 
ingly, such a feature is also present in the stress tensor of 
quantum fields on the Schwarzschild background Q. 

We found that the system of effective BW equations 
is closed to our level of approximation and leads to the 
collapsing dust star emitting a flux of energy which, at 
relatively low energies, approaches the Hawking behav- 
ior when the (apparent) horizon is being formed (let us 
note that similar features seem to appear for a quantum 
black hole [3(| as well as in the semiclassical treatment 
of collapsing shells Although we cannot determine 

a precise value for such a flux, which depends on the 
strength of the dark energy U ~ fi, consistency of the 
model constrains \i for astrophysical objects both below 
and above. With that, we were able to suggest the new 
stronger bound (|99|l for the brane tension by compar- 
ing our results with standard four-dimensional quantum 
computations of the Hawking flux for astrophysical ob- 
jects. Further, inside the star U is negative, so that each 
dust shell mostly releases energy into the next shell of 
larger radius and the whole process occurs mainly on the 
brane. This behavior then changes gradually moving to 
the exterior of the star, where U becomes positive and 
the energy lost from the core is mainly converted into 
bulk gravitational waves. 

We have also shown that the collapsing core will reach 
a minimum after a finite proper time and the collapse 
will then turn into an explosion which drives the whole 
system back to the initial state. This happens because 
the BW correction to the matter stress tensor acts as 
an "anti-evaporating" contribution which becomes big- 
ger as the energy increases. The bounce will occur after 
the formation of the apparent horizon (so that a distant 
observer presumably experiences the explosion only after 
a very long amount of time) and will not allow the for- 
mation of a global event horizon. Interestingly, from the 
Quantum Gravity side, it seems that a similar scenario 
would solve the information loss problem [38j. In fact, 
such a behavior for the core was previously obtained in 
an improved semiclassical treatment of the OS model in 
Ref. |32j , where quantum gravitational fluctuations were 
shown to have effects like those which the Weyl term 
causes in the present context. In any case, one might 
reasonably question that the bouncing ends back to the 
exact initial state. Let us then remark that matter in the 
OS model is frictionless dust, and that, in a more realistic 
case, friction would of course dissipate energy and make 
the evolution irreversible, as we discussed in Section Ivl 

The trace anomaly of four-dimensional quantum field 
theory on the Schwarzschild background has also been 
naturally interpreted as the BW correction to the trace 
of the matter stress tensor at the boundary of the core. 



Moreover, it has been shown that the back-reaction 
on the brane metric effectively annihilates the anomaly 
throughout the transition region and into the tail, com- 
patibly with the effective four-dimensional Einstein equa- 
tions, unlike semiclassical computations in which the Ein- 
stein equations are solved at the purely classical level 
(zero order in the Planck constant). Thus, if one believes 
in the holographic interpretation, it seems that the quan- 
tum anomaly would disappear to first order in the Planck 
constant when properly considering the back-reaction on 
the metric. 

Let us finally point out that all the above features were 
obtained for black holes formed by gravitational collapse, 
excluding therefore primordial black holes about which 
we have nothing to say. 
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APPENDIX A: DIVERGING EFFECTIVE MASS 
IN A TOLMAN CORE 

We shall here show that a general Tolman metric 
should have a bounce as well as an "anti-evaporating" 
phase in the BW. We shall just consider cases in which 
the space-time is globally hyperbolic and has a non- 
compact Cauchy surface (as it seems reasonable for a 
physical gravitational collapse). 

Assume that the Weyl tensor is zero. Given a null 
vector k^, from the effective Einstein equations IjlOl) we 
have 

k» k v = T M „ k v = (p cff + p cS ) (u M k^f 

This condition ensures that the singularity theorem of 
Ref. [3ij holds and the space-time will therefore reach a 
singular point at finite proper time. In particular, for the 
Tolman geometry, this means that R — > after a finite 
amount of proper time and the total effective mass in 
Eq. © with U = 0, 

M = m p + — / p 2 (R 3 )' dx = m p + AM , (A2) 
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will correspondingly diverge, as we shall now show ex- 
plicitly. 

First of all, let us note that the equation of motion for 
dust shells inside the core (0 < r < ro) yield 



R 2 R 2 



2 m p (r) + 2 AM (r, r) > 2 m p (r) , (A3) 



since AM(r, r) > 0. On considering that the flat Tolman 
solution [2J| satisfies the equation 



E?2 n2 



2m p (r) 



we then have 



i?(r,r) < R-j(t, r) 



(A4) 



(A5) 



for collapsing solutions with R < and i?x < 0. By 
re-scaling the coordinate r, one can always write |30| 



i? T (T, r) = g(r) [r (r) - t] 2/3 



(A6a) 



with g(r) as in Eq. 149(1 . and the corresponding bare mass 
is given by 

3 

(A6b) 



m p (r) = M s — 



ro 

where Ms = m p (ro). The function To(r) represents the 
(proper) time at which the shell of comoving radius r hits 
i?x = in GR and must be monotonically non-decreasing 
in r in order to avoid shell crossings (for the OS case, one 
has To = T independently of r). The BW correction to 
the effective mass can thus be estimated as 



AM = 



1 



8ttA 



R 2 R' 



■dx > 



1 



8ttA 



dx (A7) 



The last integral is not well-defined for all r > 0. In fact, 
i?x — > for t — > To(r) and AM (which is proportional 
to p 2 ) therefore diverges for r — > t ( ^" 1 (t). To see this 
more clearly, let us define T £ = tq(e) the time at which 
a shell infinitesimally close to the center (at r = e with 
< e <C ro) hits the singularity, Rt(T € ,£) — 0. For 
T > 2e, one can formally split the last integral in Eq. (|A7|) 
into two parts (we include irrelevant numerical factors in 
( ^ and finite), 



AM > C 



l ( r ) ( 



dx + C 



a _I W 



dx (A8) 



The first integration is over matter already collapsed 
into the point-like singularity (whose volume clement 
i? T R' T = 0) and the corresponding r does not represent 
a valid spatial coordinate any longer. One might try to 
regularize this integral. However, the second integration, 
which is instead over a valid range of the coordinate r, 
diverges and leads to 



AM > 



W 
Mr)-r dy 

y 



t] [3 ro (a;) + 2 x t' (x) — 3 r] 
oo . (A9) 
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FIG. 8: Error iK> II evaluated on the horizon at r = t§ s for 
M 3 = po = 1, e = 10" 4 , T = 10 and /x > (U c = 4824. The 
dashed line marks the limit of validity A < 0.5. 




FIG. 9: Error Kill evaluated on the horizon at r = t§ s for 
po = 1, e = 010 -4 , T = 10 and p. = 5000 > 4824 = /i c (M s = 
!)• 



The only way to avoid the above divergence is to have 
a negative Weyl tensor. If U ^ — p 2 /2A, the singularity 
theorem is violated and the collapse will experience a 
bounce as discussed in the text, whereas for U = —p 2 /2X 
the NLCE's imply that 



n 



(i 



XR 3 



o 2 R 2 Rl dx 



(A10) 



o 



where the dots stand for harmless terms. As we have 
shown, the above expression diverges as soon as any shell 
at r = e <C ro approaches the singularity, thus making 
the anisotropic stress tensor singular in an extended re- 
gion r > 0. If one requires that the four-dimensional 
space-time is regular everywhere and at any time, a part 
from isolated points, one should then consider U < and 
sufficiently large (with U ^ — p 2 /2X), so as to make the 
collapse bounce in the Tolman case as well. 

In light of the above analysis, we believe that the 
bouncing is a general feature of the gravitational col- 
lapse in the BW, although a numerical analysis of the 
five-dimensional equations is needed to ensure the ab- 
sence of other singularities in the bulk. 



17 



A 




. 01 



the error produced by truncating to first order in e. In 
particular, it is clear from Fig. [S] that A grows linearly 
with fj, and, from Fig. that it instead decreases with 
increasing M§. 



20 40 60 80 100 

T 

FIG. 10: Error @D evaluated on the horizon at r = r§ s for 
M s = po = 1, e = 10~ 4 , and fi = 5000 > 4824 = fi c . 

APPENDIX B: ERROR ESTIMATES 

We shall here display a few plots to clarify the behavior 
of the function A defined in Eq. (|61|l as an estimate of 



We finally show in Fig. QJJI that the approximation on 
the horizon becomes worse for increasing T, the proper 
time at which the OS core would hit the central singular- 
ity in GR. This parameter has no physical meaning for 
the bouncing core, hence can be fixed by minimizing the 
error A in the time interval of interest. Note, however, 
that we need T > 4 Ms/3 in order to have t§° > 0, and 
we also want that the horizon forms a relatively long time 
after the system begins to evolve. We therefore start this 
graph from T = 10, which is the same value we use as a 
fair optimization for the quantities plotted in Figs. 
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to dMn/dr, if it were not negligible, it would shift the 
curve in Fig. E] upwards, thus making the energy flux of 
a BW star differ more from the Hawking behavior after 
the horizon has formed. 

[60] We recall that M tx R. 

[61] In particular, for this range of masses, it was proven in 
Ref. [33] ] that only massless spin 0, 1 and 1/2 particles 
contribute. 

[62] We recall that fl% = -8ttT in GR, whereas i? M M = 

-8 7rT cff in the BW. 
[63] One would also have a divergence for r = To(r) + 

(2/3) r T~o(r), but since rd(r) > 0, the latter occurs at 

later times and we have ignored it. 



